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^3 ■ Pseudospin symmetry in nuclei is investigated by solving the Dirac equation with Woods-Saxon 

' scalar and vector radial potentials, and studying the correlation of the energy splittings of pseu- 

dospin partners with the nuclear potential parameters. The pseudospin interaction is related to a 
pseudospin-orbit term that arises in a Schroedinger-like equation for the lower component of the 
' Dirac spinor. We show that the contribution from this term to the energy splittings of pseudospin 

\^ , partners is large. The near pseudospin degeneracy results from a significant cancelation among the 

different terms in that equation, manifesting the dynamical character of this symmetry in the nu- 
cleus. We analyze the isospin dependence of the pseudospin symmetry and find that its dynamical 
character is behind the different pseudospin splittings observed in neutron and proton spectra of 
nuclei. 
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PACS numbers: 21.10.-k, 21.10.Hw, 21.60.Cs 



INTRODUCTION 



I The idea of pseudospin was introduced to explain the quasi-degeneracy in some nuclei between single-nucleon states 
r"! . with quantum numbers (n, ^, j = ^ + 1/2) and (n — 1, £ + 2, j = £ + 3/2) where n, ^, and j are the radial, the 
orbital, and the total angular momentum quantum numbers, respectively. These levels have the same "pseudo" orbital 
^ , angular momentum quantum number, £ = ^ + 1, and "pseudo" spin quantum number, s = 1/2. For example, for 
^ ' ['^'5 1/2 7 ~ l)<^3/2] oiie has £ = 1, for [np3/2, — l)/5/2] oiie has £ = 2, etc. Pseudospin symmetry is exact when 
doublets with j = £ ±s are degenerate. 

Since the proposal of the pseudospin concept, about 30 years ago [|l|, |j, the underlying mechanism of pseudospin 
symmetry in nuclei has become a topic of intense theoretical research. In a recent letter [^, we have shown the 
dynamical character of this symmetry in nuclei, arising from balance effects among the central nuclear potentials 
parameters. We have also concluded that this dynamical behavior of the symmetry is the main explanation for the 
isospin asymmetry of the nuclear pseudospin orbit interaction observed in nature. In this paper we elaborate on the 
main topics of that letter, extending our analysis for different nuclei. 

It is worth reviewing the main results of several papers [§, |[ 0, |[ ^ in which the possible underlying mechanisms 
to generate pseudospin symmetry have been discussed. 

A helicity unitary transformation of a non-relativistic single-particle Hamiltonian was considered by Blokhin et al. 
m . They showed that the transformed radial wave functions behave as 



^ ^/^^' (^\ 

^njl I ^-(£+3)^ ^^^^ W 

Moreover, they claimed this behavior to be universal and conjectured that, at small distances, the helicity transformed 
mean field of a heavy nucleus should not be much different from the untransformed one, since the behavior of R 
remains similar to R near the origin. At larger distances, however, as suggested by the above asymptotic behavior 
of the transformed wave function, the helicity transformed mean field would acquire a more diffuse surface. At the 
same time, the helicity transformed fields would get strongly nonlocal at the surface. The application of the helicity 
operator to the non-relativistic single-particle wave function maps the normal state (^, s) onto the "pseudo" state 
(£,5), while keeping all other global symmetries |4j. 

A similar kind of unitary transformation was also considered earlier to discuss the pseudospin symmetry in the non- 
relativistic harmonic oscillator It was shown that, if the ratio of the coefficients of spin-orbit and orbit-orbit 
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terms in the single-particle spherical Nilsson Hamiltonian has a particular value, pseudospin symmetry will be exact. 
Bahri et al. [|[ showed that this value was approached by relativistic mean field (RMF) estimates for that ratio. 

More recently, the subject was revived when Ginocchio j?] recognized, for the first time, the relativistic character 
of the symmetry. He noted that the pseudo-orbital angular momentum is just the orbital angular momentum of the 
lower component of the Dirac spinor. He also showed that pseudospin symmetry is an exact symmetry for the Dirac 
Hamiltonian with an attractive scalar potential, and a repulsive vector potential V, when these potentials are equal 
in magnitude: 6' + ^ = 0. It turns out that this symmetry is nothing but a SU(2) symmetry of the Dirac Hamiltonian 

In RMF models, often referred to as Quantum Hadrodynamics [l^, the nuclear saturation mechanism is explained 
by the cancelation between a large scalar and a large vector fields [|13|, [l^ , [l^, Q |l^ . Typical values for these fields 
in heavy nuclei are of the order of a few hundred MeV (with opposite signs), their sum providing a binding potential 
of about —60 MeV at the center. Therefore, it was natural for Ginocchio to regard the quasi-degenerate pseudospin 
doublets in nuclei as arising from the near equality in magnitude of the attractive scalar, and the repulsive vector, 
y, relativistic mean fields, ~ — V, in which the nucleons move. 

After the work of Ginocchio, Meng et al. [^, ^ showed that pseudospin symmetry is exact when dYj/dr = 0, where 
= S -\-V . They also related the onset of the pseudospin symmetry as a competition between the centrifugal barrier 
and the pseudospin-orbit potential. In short, they claimed that the following inequality assures pseudospin symmetry: 

where the quantum number obeys /^(/^ — 1) = i{i -\- 1) and E is the binding energy. The above inequality is, 
of course, satisfied in the particular cases E = and dT^/dr = 0. One may note that the latter condition implies 
the former if S*, V ^ at large distances, which is the case of the nucleus. However, this criterium presents several 
practical problems that are not easily overcome [|l^ , namely it cannot be applied in the range of values of r for which 
E ^ This is related to the fact that the pseudospin-orbit term has the denominator E — T,{r) and thus becomes 

infinite when E = E( r) j p, p^Sj . 

In other works [|l9|, polf 2lJ[2^|2^ the pseudospin symmetry has also been discussed in relation with the arguments 
and conjectures of refs. [^rj^^]^- Recently, a test of nuclear wave functions for pseudospin symmetry was done in [ p4| . 
In a more recent work, the structure of radial nodes occuring in pseudospin levels and a classification for the intruder 
levels (states which do not have a pseudospin partner in the hmit of pseudospin symmetry) was explained as a direct 
effect of the behavior of nodes of Dirac bound states ^ , giving support for the relativistic interpretation of nuclear 
pseudospin symmetry. However, despite these reports on the quasi-degeneracy of pseudospin doublets, the nature of 
pseudospin breaking in nuclei is not fully understood. 

In this paper we address this question by studying the splittings of the neutron single-particle energy levels for 
pseudospin partners obtained by solving the Dirac equation with vector and scalar potentials of Woods-Saxon type. 
From the behavior of these splittings with Woods-Saxon parameters, we are able to show that pseudospin symmetry 
has a dynamical character. As we have already pointed out, the present paper is a natural extension of |^], providing 
more technical details and complementary physical information. In addition, we calculate the expectation value of 
the pseudospin-orbit coupling and find that it correlates well with the observed pseudospin splittings. 

The paper is organized as follows: in Section II we review the formalism of the Dirac equation with scalar and 
vector potentials, its relation to pseudospin symmetry and identify the pseudospin-orbit term. The neutron energy 
levels and respective pseudospin splittings obtained in a relativistic calculation with Woods-Saxon mean fields are 
presented in Section HI. We consider the cases of ^^^Pb, ^^Ca and ^^Ca and draw the systematics for the variation 
of pseudospin splittings with the Woods-Saxon potential parameters. We also study the contribution to the energy 
splitting of pseudospin partners by each term in the Schroedinger-like second order differential equation for the lower 
component of the Dirac spinor, in particular the contribution of the pseudospin-orbit term. Using this systematics, 
in Section IV we explain the differences of those splittings for neutrons and protons and why this does not happen 
with spin-orbit splittings. In Section V we give a brief summary of our main conclusions. 

II. DIRAC EQUATION AND PSEUDOSPIN SYMMETRY 

The Hamiltonian of a Dirac particle of mass m in an external scalar, and vector, potentials is given by 

H = cx-p^(3{m^S)^V , (3) 

where a and p are the Dirac matrices 
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and cr are the Pauli matrices. 

The Hamiltonian (||) is invariant under SU(2) transformations for two cases: S = V and S = —V jl^, |l2[. The first 
case was recently invoked to explain the spectrum of some mesons, exhibiting small spin-orbital splittings |2^]. The 
second situation would lead to exact pseudospin symmetry in nuclei, with S = 0. However, since H plays the role of a 
binding potential, this value for S would imply that there would be no bound nucleon states 0. In relativistic mean 
field models the attractive scalar and the repulsive vector parts of the nucleon potential are similar in magnitude but 
they are not equal, and so that limiting case is never reached. 

We define /S. = V — S and denote the upper and lower components of the Dirac spinor by ^± = The 
Dirac equation = with H given by (||) and energy e = E -\- m gives rise to two coupled first-order differential 
equations for the upper and lower components: 

a.p^_ = (£;-E)^+ (5a) 
a-p^+ = (£; + 2m-A)^_. (5b) 

If one multiplies these equations by a • p one gets 

j9^^_ = -[a-pE]^+ + (E-E)a-p^+ (6) 
= -[a -pA]^- + (£; + 2m- A)a -p^- , (7) 

where the square brackets mean that the operator a • p only acts on the potential in front of it. Using Eqs. (||) one 
obtains 

= -[a-pq^+ + (£; + 2m- A)(£;-E)^_ (8) 

= -[a-pA]^_ + (£; + 2m- A)(£;-E)^+ . (9) 

Assuming that the potentials E and A are radial functions, and using again Eqs. (||), we obtain the following two 
decoupled second-order differential equations for ^±: 

= --^f -|- + -o--i:W_ + (E + 2m-A)(E-E)^_ (10) 



E — Ti\ dr r 

pH, = -^^^(-|: + i..x)*, + (E + 2.n-A)(F-E)*,, (11) 

where the primes denote derivatives with respect to r. The spinor s ^± can be factor ized into radial and angular 
parts: ^+ = \ Gi{r) ^^(^, </>) and ^_ = —Fi(r) ^~ {0, 0), where i denotes the quantum numbers of the single particle 
state. In terms of the radial functions Gi and Fi for upper and lower components, these equations become 

V'F, + (^F; + + (^ + 2m - A){E - E)F, = (12) 

^'^^ + ^ + 2m-A (^^ + ^^^) + + 2m - A)(^ - S)G. = , (13) 

where the property cr • L^f = — (1 ± was used. The cr • L terms in the equations for Fi and Gi are related to 

pseudospin-orbital and spin-orbital couplings respectively [p[ pj 

Let us look more closely at the Eq. (O) for The denominator E — in the cr • L term comes from the 



replacement of using Eq. (5a), in the term —[a • of Eq. (|8D, since 



^1 = = - — + -tr • L . (14) 

From Eq. (|5a|) , when — E = 0, then a - p^- = 0, meaning that, in Eq. ([l^), the whole term in parenthesis applied 
to ^_ gives zero. Therefore, although each of the two terms within the parenthesis, divided bv £^ — E, is infinite when 
£^ = E, their sum is finite. The analysis of these terms is better done by rewriting Eq. (g) as a Schroedinger-like 
equation: 



(15) 



4 



where an energy- and r-dependent effective mass m* = + 2m — A)/2 was defined. Applying again Eq. (|lj) we 
have 

^ + - — + -o- -L U_ +E^_ = . (16) 



2m* 2m* E — T^\ dr r 

Now we can identify the various terms contributing to the binding energy E and, in particular, assess the role of the 
term 1]' /[2m*{E — S)]cr • L/r^- in the pseudospin energy splitting. This is possible since the principal value of the 
integral of this term, after multiplying it by ^L, is finite. Furthermore, the following sum rule should be satisfied: 



J "2m*£;-S dr J " 2m* £; - E r J ~ 2m* ^ 



(17) 



where 'P' denotes the principal value of the integral. The analysis of the contribution of the various terms in Eq. (pi 
will be done in the next section. 

It is interesting to discuss how the pseudospin symmetry gets broken. The SU(2) generators of the pseudospin 
symmetry are given by [|lO|, |ll|, |l2[ 



Si 



, s, 
where 



(18) 



Si = Si — - Si, (19) 

and 5i = (7^/2 (z = 1, 2, 3). The commutator of this operator with the Hamiltonian (^) is 

[H,Si] = (^^\'^^ o) . (20) 

Thus, the pseudospin symmetry breaking can be related to the commutator If we take the non-relativistic 

limit, we recover the findings of Ref. as long as Si commutes with the non-relativistic Hamiltonian Hnr- This is 
equivalent to consider the helicity transformed spin operator, s^, to be the generator of a SU(2) symmetry of the 
non-relativistic Hamiltonian, [Hnr^Si] = 0, implying [Hnr^Si] = 0, where 

Hfir — Hfir (^-^) 

p p 

is the helicity transformed non-relativistic Hamiltonian of Ref. Requiring [S, 5^] = is equivalent to the previous 
conditions ^ = S -\- V = or dT^ / dr = 0, when S is a radial function 0, ||, |9). 



III. PSEUDOSPIN ENERGY SPLITTINGS FOR WOODS-SAXON MEAN FIELDS 

As stated above, the equivalent conditions = S -\-V = and dT^/dr = are not met in nuclei. From Ginocchio's 
findings one expects that there is a correlation between the pseudospin splitting and the S depth. Our purpose 
here is to go further and investigate whether there is also a correlation between those splittings and the shape of the 
mean-field nuclear potential, namely its radius and surface diffuseness. 

To this end, we perform a model calculation using a Lorentz structured relativistic potential of Woods-Saxon type in 
the Dirac equation. The scalar and vector components of this potential are the mean-field central nuclear potentials. 
The potentials are given by 

= l + eJ:-R)/a] ' 

where U (r) stands either for the vector or for the scalar potential. Although this potential is not a full self-consistent 
relativistic potential derived from meson fields, it is realistic enough to be applied to nuclei. Indeed, most self- 
consistent potentials have a Woods- Saxon-like shape, i.e., one can recognize in them a depth, [/q, a radius (range). 
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and a diffusivity, a. Therefore, the study of pseudospin partners sphttings as a function of these parameters is 
meaningful and reahstic enough to be apphed to most nuclei, at least qualitatively. It is known that, in certain isotope 
chains, the central depth and the surface diffuseness changes sensibly [|2^, |2^. At the same time, as the mass number 
A increases, the nuclear radius increases according to the A^^^ law, which means that it is also important to study 
the role of the parameter R in pseudospin symmetry. 

Using Woods-Saxon potentials for E and A, we solved numerically the coupled first-order Dirac equations for the 
radial fields Fi{r) and Gi{r) obtained by factoring out the angular wave functions of ^± in Eqs. dq). They read 



F~ 



G f 



1 + 







-Gt 



-{E -m-^)Ge 
{E + m- A)F^ 



(23a) 
(23b) 



where 



I j=e-l 



-ii+i) j = i+\ 



I 



(24) 
(25) 



There are altogether six parameters for E and A, namely the central depths, Eq and Aq, two radii and two 
diffuseness parameters. We observed that the pseudospin splitting is not sensitive to R and a of the A potential, and, 
accordingly, set the same radius, i?, and surface diffuseness, a, for both potentials. We first fit these parameters to the 
neutron spectra of ^^^Pb. The quality of our fitting is showed in Fig. 1, where the results of the present calculation 
(WS) are shown, together with the experimental values and those obtained using a relativistic mean field approach 
(model Gl of Ref. ||2^) for the same set of pseudospin doublets: (liii/2 , 2^9/2), (2/5/2 , 3^3/2) and (I/19/2 , 2/7/2), 
which are the 3 topmost pseudospin partner levels in ^^^Pb. 



208 



Pb neutrons 




-20 



FIG. 1: Calculated neutron single-particle energy levels of the pseudospin partners (l/i9/2, 2/7/2), (2/5/2, 3^3/2) (1^11/2, 2^9/2) 
in ^°^Pb (WS). The leftmost values are the corresponding values of the model Gl of Furnsthal et al. |29[ | and the experimental 
values are in the rightmost column. The parameters for the Woods-Saxon potential used to fit the Pb neutron energy levels 
are R — 7 fm, Ao = 650 MeV, So = —66 MeV and a — 0.6 fm. The experimental values were taken from Ref. |30| ]. 

The parameters used for the fit are given in the caption of Fig. 1. For the present purpose, our fitting is quite 
reasonable. Keeping Sq, Aq, and R fixed, we vary a in order to see how the energies of the pseudospin doublets are 
sensitive to the surface diffuseness. This dependence is shown in Fig. 2. As a increases, the splittings of the pseudospin 
doublets decrease. This could be expected because the derivative of S decreases when a increases. However, by 
further increasing a , the pseudospin doublet partners cross each other, inverting the sign of the energy splitting. This 
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inversion of pseudospin partner splittings, [+1/2 ^ ^ni-1/2 changing to ^ 1-1/2^ observed 

experimentally and was also found in ||, |8[ [l9|. It also occurs for calcium isotopes, analyzed below. 

It is important to note that, once the pseudospin doublet partners cross each other, i.e., the pseudospin splitting 
changes sign, the effect of increasing a is to drive the pseudospin doublets further apart. This systematics is consistent 
with the non-relativistic prediction of Ref. [|| as shown by the asymptotic behavior of the radial wave function (|l|). 

Pb neutrons 
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FIG. 2: pseudospin energy splittings in ^°^Pb as a function of the diffusivity for the neutron pseudospin partners of Fig. 1. The 
vertical line corresponds to the fitted Woods-Saxon parameters for ^"^^Pb given in the caption of Fig. 1. 




FIG. 3: pseudospin energy splittings in Pb as a function of the radius for the same neutron pseudospin partners of Fig. 1. 
The vertical line corresponds to the fitted Woods-Saxon parameters for ^°^Pb given in the caption of Fig. 1. 

A similar trend is seen when we let the radius R vary, fixing all other parameters. The result is shown in Fig. 3. 
When the radius increase, the energy splittings increase, which again can be understood by the dependence of dT^/dr 
with R. However, we found that the dependence of the pseudospin splitting on R is the opposite for the deepest 
energy levels, which means there are also some surface effects in this trend. This does not happen when the diffusivity 
changes, the trend being the same for all levels. Here we can see again the phenomenon of level crossing for the 
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pseudospin partners (2/5/2 , 3^3/2). For < 7 fm we could not get bound solutions for one or both liii/2, 2^9/2 
levels. This dependence on R is especially important when comparing different isotopes, since, in many cases, it is the 
radius of the mean-field potentials that changes most noticeably. Finally, we keep a, R and Aq fixed but vary Eq in 
order to study the sensitiveness of the pseudospin doublets with the depth of the central S mean field potential. The 
results are presented in Fig. 4 and the following behavior is observed: as | So | decreases, the splitting also decreases. 
This is in accordance with Ginocchio predictions for pseudospin symmetry breaking due to the finiteness of the E 
mean field. 

Pb neutrons 




_A \ \ \ \ \ L_ 

-66 -64 -62 -60 

(MeV) 



FIG. 4: Pseudospin energy splitting in ^^^Pb as a function of the depth of the S potential for the same neutron pseudospin 
partners (similar to Fig. 1). The vertical line corresponds to the fitted Woods-Saxon parameters for ^°^Pb given in the caption 
of Fig. 1. 

From Fig. 4 we see that the decreasing of \T^q\ shifts the levels (liii/2, 2^9/2) to the continuum. We also found 
that, for a deeper pseudospin doublet, an inversion of pseudospin energy splitting occurs for sufficiently low |Eo|. 



Pb neutrons 




FIG. 5: Pseudospin energy splitting in ^°^Pb as a function of R keeping a and the product |XIo|i?^ fixed for the same neutron 
pseudospin partners as in Fig. 1. 
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Actually, we have found that \T^o\ and R are not completely independent parameters. Varying and R but 
keeping a and the product |So|^^ fixed, the pseudospin splittings remain almost constant, as shown in Fig. 5. 

We observed that varying the other free parameter, Aq, does not qualitatively change the splittings. Since the A 
potential is related to the effective mass of the nucleons, the basic effect is to slightly change the nucleon separation 
energies, and especially the spin-orbital coupling, as discussed below. 

Although the Woods-Saxon potential is expected to be a good approximation for central potentials of heavy nuclei 
like lead, some calculations indicate that it might also be good enough for lighter nuclei such as calcium [^. Therefore, 
in order to show the importance of the surface diffuseness and of the radius in the dynamical pseudospin symmetry, 
we applied our model also to two calcium isotopes, where the pseudospin splitting is also small. We fitted our free 



Nucleus 


R 


a 


So 


Ao 


Level 


-E'calc 


-E'exp 




3.8 


0.8 


-78 


770 


2Si/2 
l0?3/2 


-16.0 
-15.1 


-18.2 
-15.6 




4.0 


0.8 


-64 


650 


2Si/2 
lc?3/2 


-12.6 
-12.2 


-12.4 
-12.4 



TABLE L Calculated and experimental energy values of the neutron pseudospin partners 2s 1/2 and Ids/ 2 in °Ca and ^Ca. 
The Woods-Saxon potential parameters used to fit the energy levels of these nuclei are also displayed. The radius, R, and the 
diffuseness, a, are in fm and the energies in MeV. The experimental values are from Ref. |30|]. 



parameters to the experimental values of the energies of the topmost neutron levels of ^^Ca and ^^Ca. The results for 
the pseudospin partners (2si/2, l<^3/2) are presented in Table 1 and compared with the experimental values. 

The difference between the ^^Ca neutron potential with respect to the ^^Ca neutron potential lies basically in a 
smaller modulus of the central depth and in a larger radius. We also see that the magnitude of pseudospin energy 
splitting decreases from 0.9 MeV to 0.4 MeV as the number of nucleons A increases. Since the levels 2si/2 and 
lds/2 are inverted in ^^Ca, the effect of the change in R predicted by the systematics is to increase the difference 
£^2si/2 ~ ^1(^3/2' since this value is negative, to decrease its magnitude, favoring the symmetry. 

We are going now to examine more closely the effect of the term Ti'/[{E — E)(2m*)]cr • L/r^- in Eq. ( p!^ ) on 
pseudospin energy splittings. If we multiply that equation by ^L, integrate, and then divide the result by J ^L^-d^r, 
we get 



2m' 



(^pso) + (^d) + (S)=E 



(26) 



where 



2m* 



(^pso) 







\ dr J 




J 





dr Ff 





rOO -1 y(/ 

F rdr im:^2 



f 

Jo 



2m* E-H 



dr 



dr 



1 ^ dF, 



2m*£;-S ^ dr 



dr F^ 



r^ dr E F^ 



r^ dr F^ 



(27) 



(28) 



(29) 



(30) 



These terms can be identified, respectively, as a kinetic term, a pseudospin-orbit term, a potential term related to 
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what is sometimes called Darwin term and the mean value of the E potential with respect to the lower component, 
We computed numerically the principal values, checking that the sum rule (|l^) was satisfied. 

Pb neutrons 




FIG. 6: Contribution of the terms ( [2 7] ) - (|30|) to the energy splittings of the pseudospin partners of Fig. 1. The strong 
cancelation among the terms produces a small AE. 

The results are shown in Figs. 6, 7 and 8. In Fig. 6 is plotted the difference of the terms in Eq. ( p6|) between 
each member of the pseudospin partners we have been considering. One sees that the contribution for the pseudospin 
energy splittings of the Vpso term is larger than the splitting itself and it is at least of the same order of magnitude 
and has the opposite sign of the kinetic and (S) terms. The relative contribution of each term depends on the 
particular pseudospin doublet we are considering. This can be seen clearly in figures 7 and 8, where the differences 
of the terms in ( p6| ) for the pseudospin partners (l^ii/2, 2^9/2) and (2/5/2, 3^3/2), respectively, are plotted. For 
states with higher angular momentum the magnitude of (Fpso) contribution for the energy splitting is higher, since 
it is roughly proportional to 2i + 1, but there is a significant cancelation with the kinetic and (E) terms. A similar 
finding was already reported in [^, where it was pointed out that pseudospin-orbit coupling cannot be treated as 
a perturbative quantity. One the other hand, these figures show a clear correlation between the pseudospin-orbital 
term and the pseudospin energy splitting when the diffusivity and the depth of S potential are varied. This feature 
does not depend on the particular pseudospin doublet under consideration. 

The pseudospin-orbital term (Fpso) comes from the coupling between the upper and the lower components. This 
term will vanish only when the upper component is zero (assuming a non-zero lower component) or when = 0, 
showing again that exact pseudospin symmetry is possible only when there are no bound states. This is not the case 
for the usual spin-orbit interaction generated by the term in the square brackets in Eq. Here, this term will vanish 
when the lower component is zero, which corresponds to the non-relativistic limit. 

From the results presented so far, we can conclude that the degeneracy of pseudospin levels is very much dependent 
on the shape of the nuclear mean field potential: the actual choice of the potential parameters, fitted to describe a 
nucleus, determines whether or not there is pseudospin symmetry. Furthermore, it is the strong cancelation presented 
in Figs. 6-8 that produces the quasi-degeneracy of pseudospin levels. The conclusion is that pseudospin symmetry is 
a dynamical symmetry in nuclei, in the sense of Arima's definition of dynamical symmetry [|2^: (i) a symmetry of the 
Hamiltonian which is not geometrical in nature; or (ii) a ordered breaking symmetry from dynamical reasons. This 
dynamical nature of pseudospin symmetry is further stressed by the fact that pseudospin levels can cross each other. 
We observe also that the realization of pseudospin symmetry depends on the doublet in consideration, although the 
systematics is the same for all pseudospin partners. 



IV. ISOSPIN ASYMMETRY OF PSEUDOSPIN SPLITTINGS 



We apply now this systematics to some nuclei studied in the literature, looking at the differences between neutron 
and proton spectra as far as pseudospin symmetry is concerned. 
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FIG. 7: Contribution of the terms - ( ^q|) to the energy spUtting of the pseudospin partner (liii/2, 2^9/2) in ^°^Pb as a 
function of the diffuseness. The vertical line corresponds to the fitted Woods-Saxon parameters for ^°^Pb given in the caption 
of Fig. 1. 
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FIG. 8: Contribution of the terms ( |27| ) - ( |30| ) to the energy spUtting of the pseudospin partner (2/5/2, 8^3/2) in Pb as 
a function of So- The vertical line corresponds to the fitted Woods-Saxon parameters for ^°^Pb given in the caption of Fig. 1. 



Recently, the pseudospin symmetry in Zr and Sn isotopes was investigated as a function of the number of nucleons 
The form of the E potential as a function of the radial distance for such nuclei, as A increases, is given in Ref. [ p8[ , 
starting with ^^^Sn and going up to ^^^Sn. This information on E for Sn allows us to analyze the behavior for neutron 
pseudospin doublets going from A = 100 to A = 170. As A increases, the modulus of the central depth, |Eo|, 
decreases and the surface diffuseness increases, both effects favoring the pseudospin symmetry, as shown in Figs. 2 
and 3. However, as expected, the radius R increases with A, which can (at least partially) offset the previous changes. 
This is what is observed in Ref. Q : the energy splittings of deep neutron pseudospin partners decrease but those near 
the surface do not change sensibly as A increases. The splitting for the partner (I/19/2, 2/7/2) increases from a small 
negative value to a positive value, i.e., the levels cross each other. In |^ it is also observed that pseudospin splittings 
are smaller for levels close to the nucleus surface, even crossing each other. This phenomenon can also be observed in 
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^^^Pb, as one can see from the experimental energy values of the pseudospin partners (2/5/2, 3^3/2) (l^ii/2, 2^9/2) 
in Fig. 1. We do not expect our ansatz for the nuclear central potentials to explain the behavior of levels very near 
the nucleus surface. Such fine details are outside the scope of our work, which simply intends to establish a broad 
systematics for pseudospin splittings in nuclei. 

For Zr isotopes, studied in refs. [^, ^, although the authors did not provide details on the potentials, again one 
sees that the pseudospin splitting in general decreases with A, with a somewhat irregular behavior for the partner 
(1^7/2, 2(i5/2)- It is interesting to remark that the splitting for the doublet (2^3/2,351/2) decreases to the point of 
level crossing. Therefore, we can conclude that, as A increases, there is a general trend for the pseudospin doublet 
partners to approach each other, towards realization of the symmetry. 

The proton spectra of Sn isotopes was also studied in Ref. [^. We observe that the central neutron depth, Hq, 
varies from about —65 MeV to —54 MeV and the surface diffuseness increases whereas the proton central depth, Sq, 
varies from —50 MeV to —60 MeV and the surface diffuseness slightly increases. The reason lies on the p meson 
interaction, which for a neutron rich nuclei is repulsive for the neutrons and attractive for the protons, as explained 
by mean field model calculations [|l3|. The inclusion of this interaction, which is important in asymmetric nuclei, 
changes the vector part of the potential V to 

V = V^^Vp = V^±^po, (31) 

with + and — signs for protons and neutrons respectively; po is the time component of the p field, which is proportional 
to the number of protons minus the number of neutrons, and K; comes from the vector-isoscalar co meson. The increase 
of the parameters R and a with A can be traced back to the known nuclear radius A^^^ dependence and to the excess 
of neutrons on the surface (neutron skin effect). Hence, we may conclude that, for a given nucleus, the parameters 
So, a and R for protons and neutrons are different. Then, an isospin asymmetry in the pseudospin interaction is 
expected to take place in agreement with the analysis of the mean field nuclear parameters for protons and neutrons 
done in Q. In particular, since po is negative for heavy nuclei, the vector potential ( [3l| ) is bigger (and, thus, |S| is 
smaller) for neutrons than for protons. From our previous analysis, for a neutron rich nucleus, pseudospin symmetry 
for the neutron spectrum is favored, in agreement with the results presented by Lalazissis et al [ p!9| . 

The systematics discussed here seems to be quite general and a comment on how it affects the spin-orbit splittings is 
pertinent. Spin-orbit splittings are much less changed by variations of the central depth and surface diffuseness of the 
nuclear potential than pseudospin-orbit splittings. The reason is that to change substantially spin-orbit splittings one 
needs a significant relativistic content, i.e., a large lower component in the Dirac spinor |3^. The spin-orbit splitting 



is completely correlated with the nucleon effective mass [|3^. From Eq. (13), we see that it will depend strongly on A. 
This potential itself carries a quite large scale when compared with S (around a factor ten). Therefore, the p meson 
potential, Vp, generally about ten percent of E, becomes irrelevant compared to A. This justifies why the spin-orbit 
interaction is roughly isospin symmetric, being almost the same for neutrons and for protons j2^, |3^ . 

On the contrary, the pseudospin splitting needs less relativistic content, i.e., a smaller lower component in order to 
change. This may explain why the non-relativistic analysis of Ref. [|[ of the origin of the pseudospin symmetry, and 
manifested by Eq. (0), works well to explain the small pseudospin splitting. The pseudospin orbit interaction depends 
on E, as one can see from Eq. (|l|). This means that now Vp cannot be neglected, since its values are comparable 
to those of E, causing sensibly different variations of the central depth for neutrons and protons, for which Vp has 
different signs. This is one of the mechanisms generating the isospin asymmetry in the pseudospin-orbit interaction, 
as we have pointed out before. This can be seen in Fig. 9, where the energy splittings for the pseudospin partners 
(I/19/2, 2/7/2), (2/5/2, 3^3/2) and for spin-orbit pairs (2/5/2, 2/7/2), (1/^9/2, 1^11/2) are plotted. Here we have just 
fixed the radius and diffusivity, using the values that fitted the neutron spectrum of ^^^Pb, and let both Eq and Aq 
change by the same amount, thereby simulating the increase of the vector potential V that affects E and A in the 
same way. It is clear that the splittings change much less for the spin-orbit doublets than for the pseudospin doublets. 
The magnitude of the change depends on the specific level one considers, but we checked that this kind of behavior 
occurs for deep levels. Of course, one should also consider the effect of the Coulomb potential for proton spectra, 
which has the opposite sign of Vp for protons. However, we expect the effect of Vp to be dominant for neutron rich 
nuclei. 



V. CONCLUSIONS 



We have studied, through a mean-field model calculation with Woods-Saxon potentials, the role played by the 
surface diffuseness, the radius, and the central depth of the sum of the standard vector and scalar nuclear potentials 
in the energy splittings of pseudospin partners, both for lead and calcium isotopes. This study allowed us to draw 
a systematics for the behavior of the pseudospin splitting as a, and Eq vary. Such behavior was confirmed when 



12 




—i \ \ \ \ \ \ \ L 

-66 -64 -62 -60 -58 

E„ (MeV) 



FIG. 9: Energy splitting in ^°^Pb as a function of the depth of the H potential for the pseudospin partners (1/19/2, 2/7/2), 
(2/5/25 3^3/2) (hollow symbols) and spin-orbit pairs (2/5/2, 2/7/2), (1/19/2, l/iii/2) (filled symbols). The depth Ao varies from 
639 MeV to 649 MeV so that only the vector potential changes. The vertical line corresponds to the fitted Woods-Saxon 
parameters for ^°^Pb given in the caption of Fig. 1. 

we applied the systematics to the Sn nuclei isotope chain, predicting the observed general trend towards degeneracy 
of pseudospin doublets. Based on our findings and in general features of the neutron and proton potentials, we 
explained the observed isospin asymmetry of pseudospin splittings. We were able to identify a pseudospin-orbital 
term and found that its contribution to the pseudospin energy splittings is large. The near degeneracy of pseudospin 
levels is obtained from a significant cancelation among the different terms in the Schroedinger-like equation for the 
lower component of the Dirac spinor. We also showed a clear correlation between the pseudospin-orbital term and 
the pseudospin energy splitting when the diffusivity and the depth of S potential are varied. These findings led us to 
argue that pseudospin symmetry in nuclei is dynamic. 

We believe that the systematics observed in our model calculation and the non-perturbative nature of the 
pseudospin-orbit term, which is larger than the energy splitting itself and needs to be canceled in order to pro- 
duce the pseudospin quasi-degeneracy may help to understand how the pseudospin symmetry is dynamically broken 
in nuclei. 
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